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Abstract 

At the Mott transition of a quark/meson plasma, mesons become un- 
bound. This leads to the effect of critical scattering, which is studied 
by investigating photon pair production due to the process qq 77. 
This process can proceed either via direct annihilation or via the for- 
mation of a mesonic resonance. It is shown that the latter channel 
leads to an enhancement of photon pairs with invariant mass equal to 
the thermal pion mass. The size of this effect measures the time the 
temperature stays near the Mott temperature during the evolution. 
It is particularly pronounced for a first order phase transition. The 
7]-0 —J. decay gives a strong background contribution and may make 
the observation of critical scattering in two-photon spectra of present 
day heavy ion collision experiments difficult. 
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I. INTRODUCTION 



In ultrarelativistic heavy ion collision experiments at CERN SPS and BNL AGS 
hadronic matter is produced under extreme conditions of temperature and density. 
These conditions might be sufficient for the transition to a new phase of matter, 
the quark gluon plasma ||T],3. Up to now, no unambiguous signature of this hy- 
pothetical phase has been found, and it appears that a rather complex analysis of 
observables has to be performed. The inputs for such a numerical analysis should 
be provided from effective approaches to QCD in the nonperturbative low energy 
sector where chiral symmetry breaking and confinement of QCD degrees of freedom 
are the prominent features. 

Although plagued with nonrenormalizability and the absence of confinement, 
the Nambu-Jona-Lasinio (NJL) model has been developed recently beyond the 
applications for calculation of hadronic mass spectra towards a kinetic description 
of hadronic matter at the quark-hadron phase transition . A particularly inter- 
esting problem within this context is the description of hadrons as bound states of 
quarks and antiquarks and their dissolution at finite temperature and density due 
to thermal excitation or to compression, i.e. the Mott effect 0. At temperatures 
below the Mott temperature Tm, the model contains massive constituent quarks 
and bound state mesons as degrees of freedom. For temperatures above Tm, these 
bound states dissolve and become resonant (nonperturbative) correlations. Whereas 
the appearance of free quarks at T < Tm is an artifact of the NJL model, the situ- 
ation for T > Tm is fairly realistic, except for the absence of gluons. We thus apply 
the model for the description of mesonic correlations at T > Tm- It has been argued 
elsewhere that the Mott transition leads to the effect of critical scattering and 

may have experimental consequences. In this work, we investigate the influence of 
such correlations in a quark plasma on the two-photon spectrum which is one of 



the observables measured in nucleus- nucleus collisions, see e.g. |Tl|, and which has 
been previously studied in Refs. To this end, we present a model study of 

the process qq —>■ 77 employing the two flavor NJL model. In the s-channel, this 
process proceeds via the formation of a virtual meson. At the Mott temperature, 
this meson becomes unbound, which leads to a divergent scattering length. We es- 
timate the magnitude of this critical scattering effect and its sensitivity to the order 
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of the phase transition by comparing the spectra containing resonance contributions 
with those, which only contain contributions of direct annihilation. However, we 
give only a crude estimate of the contributions from 7r° decay into photons after 
freezeout . 

Since presently no consistent transport theory modeling the chiral phase transi- 
tion is available, we apply a hydrodynamical scenario in order to obtain the total 



photon spectra. In this approach, we use a model equation of state |ll^ , p!5| which in- 
terpolates smoothly between a first and a second order phase transition. The results 
of our investigation show that the yield of photon pairs with invariant mass m7r(Tjv/) 
depends strongly on the time the temperature stays near the critical temperature. 
This paper is structured as follows: The calculation of photon production cross 



sections including resonance contributions is performed in Section 0. In Section |T] 



we apply the Bjorken hydrodynamical expansion scenario |TH] in order to calculate 



the invariant mass spectrum of photon pairs from a hadronizing quark plasma. In 



Section we summarize and present the conclusions. 



II. PHOTON PRODUCTION CROSS SECTIONS 

A. Review of Formalism 

For the calculation of photon production cross sections we employ the Nambu- 
Jona-Lasinio (NJL) model in its SUf{2) version The Lagrangian for this model 
reads 

C = ^{t^ - mog)i; + G [ii'^Y + iiji^^rtpy] . (2.1) 

Here G is a coupling constant of dimension MeV~^ and r the Pauli matrices in flavor 
space. The quark wave function ip carries flavor and color degrees of freedom, which 
are implicitly summed in Eq. The symbol rrioq denotes the current quark mass, 

which explicitly breaks chiral invariance. 

The masses of the constituent quarks are computed in leading order in an expan- 
sion in the inverse number of colors, 1/Nc |]T^. For the constituent quark masses. 



this corresponds to a Hartree approximation. One obtains the gap equation for the 
physical masses 
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niq = TTioq + 4iGNctT^S{x, x) = rrioq -^mgA{mg) , (2.2) 



where the function A(mq) is given in the imaginary time Matsubara formahsm |]TB| 
as 

|p1<a 



Here we have abbreviated E = ^fp' + m^, and u;„ = (2n + l)7r//5 is the fermionic 
Matsubara frequency. The symbol (3 = 1/T denotes the inverse temperature. The 
function A{mq) is given exphcitly by [|l^ 

A(m„) = -4 /^rfp^tanh(/5E/2) , (2.4) 
JO E 

where we have used an 0(3) cutoff to make the integral finite. 

Mesonic correlations in the pseudoscalar (pion) and scalar (sigma) channels of 
the qq interaction are treated in the random phase approximation [^]. This leads to 
the form 

2(7 

P^(A;o, k) = ^ (2.5) 

^ ^ 1 - 4Gn^(A;o, k) ^ ' 

for the pion propagator and 

2G 

VAko,k)= , (2.6) 

1 -AGIl^{ko,k) 

for the a propagator. In Eqs. ( ^.5] ) and (^.61), 11^ and 11'^ denote the irreducible pseu- 
doscalar and scalar polarization functions, respectively. They are given explicitly by 



n^(A;o, k) = [2A{mq) - {k^ - e)Bo{ko, k)\ (2.7) 



n^(A;o, k) = [2A{mq) + {Am'q - kf, + A;")5o(A;o, k)\ , (2.8) 

where the function BQ{kQ, k) is defined as the analytical continuation of 

|p1<a 
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to real values of iz/„ 



ko. In Eq. ( p.9|) , we have used the notation E' 



{p — ky + m'^. In the following, we will need the function Bq for the special 



case k = 0, where it can be expressed as [19 



(2.10) 



+ in 



\ 



4m^ 
1 tanh 

Kg 



'Pko\ 



Q{ko - 2m,) 2.,/A2 + m^- k 



One can immediately see from this form, that -Bo(^05 0) is a real function for k^ < 2m, 
and a complex function for ko > 2mq. 

The masses of the vr and a mesons are obtained via the dispersion relations [H] 



1 - 4Gn^(m^, 0) = 
1 -4G'ff(m^,0) = . 



(2.11) 
(2.12) 



Note that these are in general complex equations and the solution may contain a 
nonzero width: 



TTlr 



a 



TT, cr 



(2.13) 



In this case, the function Bolk^, 0) has to be analytically continued to complex 
values of fco ID- However, Eq. ( p.ll|) at low temperatures has real (bound state) 
solutions for m^ with mj^{T) < 2mg(T). At the Mott temperature Tm-, one obtains 
f^iriTM) = 2mq(TM)- FoY higher temperatures, the pion becomes a resonant state. 
Further details about this transition can be found in Ref. The sigma meson 
differs from the pion in that it is a resonance for all temperatures. In Fig. ^ we 
show the result of a numerical solution of Eqs. ( p.2|) , ( P . 1 1 ) and ( 2.12 ), using the 
parameter set mo, = 5.0 MeV, A = 653.3 MeV and GA^ = 2.105, as a function of 
temperature. At low temperatures, the pion mass, which is indicated by the solid 
line, lies below the mass of its constituents, which is given by the dashed line. At 
T = Tm, one obtains mT,(TM) = 2mg(TM), and the pion moves from being a bound 
state to a resonant state. For our parameter set, this happens at Tm = 205 MeV. 
The mass of the sigma meson, which is indicated by the dotted line in Fig. |l], is larger 
than 2mg in the whole temperature range. The sigma meson is thus a resonant state 
at all temperatures. 
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B. Photon Production Cross Sections 



The Feynman diagrams, which contribute to the process qq —> 77 in lowest order 
in 1/Nc are shown in Fig. Beside the direct terms, i. e. those which do not contain 
resonance contributions and which are already present in quantum electrodynamics 
I P^Of , there are four diagrams, in which the incoming quarks first form a virtual 
TT^ or a state, which afterwards decays into two photons. We write the transition 
amplitude for the t-channel direct term shown in Fig. ^ as 

- iMt = Sc,,c2V2iP2)i-iqfe^2) , _u _ . (2.14) 

In this equation, pi and p2 denote the four-momenta of the incoming quarks, ki 
and ^2 the four-momenta of the outgoing photons and ei and €2 their respective 
polarization vectors. The electron charge magnitude is denoted by e, g/ is the quark 
charge factor; one has q^ = +2/3 for up quarks and qa = —1/3 for down quarks. The 
Kronecker symbol 5ci,c2 guarantees that only color singlet states can annihilate into 
two photons. The transition amplitude for the w-channel direct graph is evaluated 
analogously to be 

1 

-iMu = Sc,,c2V2{P2){-'iqfeh)j u (-%e/2)M(pi) • (2.15) 

Multiple rescattering in the initial state leads to the formation of a strong correlation 
in the s-channel as is shown diagrammatically in Fig. ^ where the correlation is 
described by a meson propagator. The transition amplitude corresponding to sigma 
propagation in the s-channel can be written as 

-iMs,a = 5c„cAp2)u{j)i) [2l?,(v^,0)] [-ia{s){k^kP ~ g^^k^k^)] t^^e^x . (2.16) 

In this expression, we have restricted ourselves to the kinematical situation that the 
center of mass system rests with respect to the medium. The factor —ia{s){k^k2 — 
g^^kik2) stems from the triangle part of the diagram in Fig. The explicit form 
of a{s) is derived in Appendix 0. Here we only note that a{s) contains a factor 
(1 — Arrig/s), which leads to a suppression of this diagram near threshold. The 
transition amplitude corresponding to the propagation of a pionic mode is given by 

- iMs,n = 5ci,c2viP2)iil5Tff)u{pi) 0)^(^/3, 0) -ib{s)e^'" ^^k^k^ 61^62^ , (2.17) 
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where the triangle contributions are contained in the factor —ih{s)e^^ y^^k^kf^. Here, 
e denotes the totally antisymmetric tensor. As was the case for a(s), an explicit 
expression for h{s) is given in Appendix 0. Compared to Eq. ( p.l6|) , this expression 
also contains an additional isospin factor r|j, which is equal to +1 for up quarks, 
and —1 for down quarks in the initial state. 
The integrated cross section is given by 



where 



<^,.--n^r)= ^,^^(^_, ,) / dt\M\^ , (2.18) 



1 



l-^P = ^ E IM + Mu + Ms,. + MsA^ (2.19) 



c s.c 



is the spin and color averaged square of the transition amplitude and the integration 
limits are given by 



s 1 I 

- - = tmin + -^s{s - Ami) ■ (2-20) 



Explicit expressions for and its integral over t can be found in Appendix 0. 

At the Mott temperature, the pion appearing in the s-channel graph of Fig. ^(b) 
switches from a bound to a resonant state. It is a well known fact from scattering 
theory, that a bound state at threshold leads to a divergent scattering length. In the 
present calculation, this can be shown explicitly by noticing that the pion propagator 

contains a singularity at T = Tm- The leading behavior of can be obtained 
from the observation that -Bo(v^, 0) can be expanded in a power series in {Amg — sY^'^ 

BoiV^, 6)=A + B{Aml - s)i/2 + c{Aml - s) + 0((4m2 - s)'/^) . (2.21) 

The pion propagator can thus near threshold be approximated by 

2G 

l^.(v^,0) 



1 + ^ (A{mg) - 2mM) - | sB{Aml - s)i/2 -{A + sC)(4m2 - s] 



(2.22) 

with ^ = {NcG)/tt^. Since one has itIt, = 2m^ at T = Tm, the first square bracket 
in Eq. (|2.22| ) has to vanish at this temperature and one obtains 
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I)^(v^,0) ~ (s-4mJ)-i/2 _ ^2.23) 

Via Eqs. (P12|) and ( p.l8|) , one obtains 

{s - 4m5)-3/2 (2.24) 



near threshold. The pion propagator diagram thus leads to a strong enhancement 
of the production of soft photons at the Mott transition when compared with the 
direct contributions. Due to the temperature dependence of the pion mass (see 
Fig. |l|) the peak in the two-photon invariant mass spectrum at y/s = is shifted 
upwards to a value of m-^iTM) = 173 MeV. Note that the effect of a rising pion mass 
at high temperatures is also visible in other effective models [pT| , p2| and in lattice 



calculations . Moreover, the 7r° — 77 decay rate is strongly enhanced at the Mott 
transition which is an effect akin to the phenomenon of critical opalescence known 
from solid state physics . Numerically, one finds that the resonance contributions 
are most pronounced near the Mott temperature, while their importance is strongly 
reduced compared to the direct terms at higher temperatures. A similar behaviour 



has been found previously in Refs. |24,25 



III. HYDRODYNAMICAL EXPANSION 

The cross sections calculated in Sec. || are not measureable directly in heavy ion 
experiments. Instead, one has to integrate the production rates over the evolution 
of the system. A crude model for the expansion and cooling of the plasma has 
been given by Bjorken [1^. Within this model, we will use different ansdtze for the 



equation of state and investigate their influence on the physical spectra. 

A. Production Rates per Space Time Element 

The production rate of photon pairs with invariant mass M per unit space-time 
volume is given by 

= (2^(2^"-"-^"^ 
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where frix) = l/(exp(a;) + 1) is the Fermi distribution function and a^^ the sum 
of the cross sections for uu 77 and dd —>■ 77. As in Section 0, we neglect effects 
that stem from a relative motion of the quark pair with respect to the rest frame of 
the medium. In principle, Eq. should be supplemented by Bose enhancement 
factors, which account for the stimulated emission of photons in the medium. Since 
the form of such enhancement factors is unknown and can only be determined in 
a nonequilibrium theory, we discard any such factors |l^JT3| , assuming that the 
photon density stays low due to the short lifetime of the plasma and thus prohibiting 
equilibration. Note that the inclusion of an enhancement factor, even at the finally 
equilibrated level, only alters our results quantitatively by at most a factor of 10. 
Our qualitative picture remains however unchanged. The fluid velocity, which has 
to be supplied by the expansion model, is denoted by u^. As it turns out, the final 
result does not depend on u^, so one can evaluate under the assumption that 
the fluid is at rest. The final result for the production rates reads JT^pB 



dN,, (2N^YM^^M^-Amla,,{M\T) ^ l + exp(-/3E_) 

— aJ^ J F[pj^) log- 



(3.2) 

with 



(M2-2m^)^± JM2(M2-4m2)(E2-m2) . (3.3) 



An expansion of Eq. (|3.2|) near threshold shows that the singular part of the pro- 
duction rate behaves as 

At T 7^ Tm, where a^^ displays only kinematical singularities, the production rate 
thus vanishes at threshold. At the Mott temperature, however, one obtains from 
Eq. ( p.24|) , that the production rate diverges as 

dN^ 



'77 



(M2-W)-V2 (3.5) 



d^xdM 

near threshold. Note that this is an integrable singularity, i. e. although the pro- 
duction rate contains a divergence, the total number of photons produced stays 
finite. 
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B. Expansion Model 



As a simple model for the hydrodynamical expansion, we employ Bjorkens one 
dimensional model . The essence of this model is contained in the two equations 
for the fluid velocity 

1 



= -(t,0,0,z) (3.6) 
r 



and the entropy density s 

sir) To 



. ^ , (3.7) 

s{to) t 



where r = yW—z^ is the proper time of the fluid element. Note that all thermo- 
dynamic quantities depend on r only. In Eq. (pI7|), tq denotes the formation time 
of the plasma, which we treat as a free parameter in the following. To compute 
the temperature distribution from Eq. ( p. 71) , one has to make assumptions about 
the equation of state of the plasma. In the following, we will compare two different 
scenarios: 

(i) First order phase transition: The entropy density changes discontinously at 



T = Tm- In this case, the equation of state reads |114| , p!2 

do if T > Tm 

2 I <^ 



QO^ fdQ + (1 - f)dH if T = Tm . (3.8) 
dn if T < Tm 

The parameter / in Eq. ( p.8|) is the volume fraction of the system, which is in the 
quark-gluon phase, (1 — /) the volume fraction, which is in the hadron phase. The 
symbols di denote the effective numbers of degrees of freedom; one has dq = 37, 
dn = 3. From Eq. ( p. 81) , one obtains the following behavior of temperature: 

1. The expansion starts at r = tq with an initial temperature Tq. Afterwards, 
the temperature drops according to 

T = ToM , (3.9) 



until it reaches the Mott temperature at time 

' T 
■■Tm 



n = To{^)' , (3.10) 



at which the mixed phase begins. 
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2. The second stage of the expansion is isothermal with T = Tm- During this 
stage, the relative volume of the plasma evolves according to 

f{r) = (dQ^ - dn) . (3.11) 



This stage ends at time 



dq ~ dn \ T 



r2 = ^n , (3.12) 
dn 

at which f{T2) = and the system enters the purely hadronic phase. 
3. During the third stage, the temperature of the system drops as 

T = Tm i^-fj . (3.13) 

This stage lasts until breakup. 

The solid curve of Fig. ^ gives a specific example for the parameters Tq = 250 MeV 
and Tq = 1 fm/c. The temperature stays constant from ti = 1.8 fm/c to r2 = 
22.4 fm/c. 

The physical photon spectra have to be obtained by a space-time integration 
of the production rates. In the Bjorken model, the integration over the transverse 
space directions simply contributes a factor vri?^, where Ra is the radius of the 
colliding nuclei. Since the thermodynamic quantities only depend on r, it is useful 
to switch to the integration variables r and 

1 , t + z , , 

The final result is 

which gives the number of photons produced in the central rapidity interval of the 
collision. Note that Eq. (|3.15|) contains only contributions from the plasma phase. 
Although the NJL model would allow for qq annihilation in the hadronic phase, we 
discard these unphysical contributions by restricting the integration limits. 

(ii) Second order phase transition: For a second order phase transition, the jump 
in Eq. ( |3.8| ) is smeared out over an interval AT according to ||15[ 
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90 



-T 



do - dn \ AT 



(3.16) 



In this scheme, the time evolution of the temperature has to be calculated numeri- 
cally using Eq. ( p.7| ). The result is a smooth curve, which for AT approaches 
that of the first order transition, as can be seen from the dashed (AT = 10 MeV) 
and dot-dashed (AT = 1 MeV) curves in Fig. |. However, this scenario displays the 
qualitative difference to the first order scenario, that it does not develop a mixed 
phase. Note that for To — Tm ^ AT, the phase transition is reached at time 

\TmJ dg + dH 

For the parameter set of Fig. |^, this yields ti = 3.35 fm/c. This is very small 
compared to the end time of the mixed phase in the first order scenario. 
In contrast to Eq. ( |3.15| ), the photon spectrum is now given by 

-nRirdrr4^ . (3.18) 



dMdr] ^ Jro d^xdM 

As Eq. (|3.15| ), the integration is cut off at the phase transition and thus contains 
only contributions of a pure quark phase. 



C. Photon Spectra 



The physical photon spectra at central rapidity for the case of a first order 
phase transition are shown in Fig. ^ as a function of the invariant pair mass. The 
dashed line gives the photon spectra without resonance contributions. Since the 
system stays at temperature T = Tm for a finite time, the spectra develop a cusp 
at M = 2mg(T/v,/). The rise in the photon yield at larger pair masses reflects the 
contribution of the mixed phase. The photon spectra containing resonance contri- 
butions is given by the solid line in Fig. ^. Here the influence of the mixed phase is 
more dramatic: Due to the second term of Eq. (|3.15|) , the threshold singularity in 
Eq. (|3.5|) immediately carries over to the physical spectra and one thus observes a 
great enhancement of photon pairs at the invariant mass M = 2mq{TM) = itl-^ITm)- 
For our parameter set, this occurs at M = 173 MeV. Note that this peak is in princi- 
ple separated from the hadronic contributions due to tt^ decay (M = 135 MeV) and 
due to 7r+7r~ annihilation (M > 279 MeV). However, we suspect that this medium 
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effect might hardly be visible in experiment since it is not clear if the enhancement 
of the pion mass is sufficient to separate the critical scattering peak from the huge 
peak in the photon spectra stemming from the decay of "cold pions" originating after 
hadronization. A quantitative estimate of the number of photon pairs originating 
from 7r° decay can be made from the relation |T^ 

-^^ = 1.5roT^ = 3.Um-' . (3.19) 

The total number of photon pairs contained in the solid curve of Fig. ^ is 

1 ^^-^^ -4.0X lO-^fm-2 , (3.20) 



which is of the order of a factor a"^^ less than the value given in Eq. (|3.19|) . The 
number of photon pairs contained in the Mott peak is 



The Mott peak thus contains roughly 1% of the total number of photons stemming 
from qq annihilation. 

In the case of a second order phase transition, the system does not stay at a 
fixed temperature for a finite time. Thus the Mott peak visible in Fig. ^ should 
disappear if the order of the phase transition is weakened. Figure |] shows the effect 
of the parameter AT on the photon spectra. Here the case of the first order phase 
transition is shown as the solid line. The other lines show the spectra in the case of a 
second order phase transition for various values of AT. One recognizes that one still 
obtains a strongly enhanced photon yield at M = m^(TM) for very small values of 
AT. For larger values of AT, however, this peak gets flattened. At AT = 20 MeV, 
only a weak cusp is visible in the two-photon spectra from a quark-meson plasma. 



IV. SUMMARY AND CONCLUSIONS 



We have investigated the influence of resonance contributions on the process 
qq — > 77 in a quark meson plasma at temperatures larger than the Mott tem- 
perature. We have employed a chiral quark model as a microscopic approach to 
a consistent description of the chiral phase transition and the Mott effect of the 
pion as a qq bound state. We have found that at the pion Mott temperature Tm, 
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which corresponds to the chiral transition temperature present for = 0, 
there occurs a narrow peak in the two photon spectrum at the invariant mass 
M = niT.iTM) ~ 173 MeV. Note that lattice resuhs give = Tdeconfinement within 
the present accuracy It seems thus natural to identify the Mott point with the 
deconfinement temperature. 

In order to estimate the role of this critical scattering effect on the photon spectra 
emitted from an expanding quark plasma, as produced e. g. in heavy ion collisions, 
we have used the hydrodynamical Bjorken scenario. The expansion dynamics fixes 
the temperature near T = Tm for a finite time interval, depending on the parameter 



AT of Eq. (|3.16|) , which is a measure for the smoothness of the transition. In 
the case of a first order transition this leads to a divergence in the photon spectra 
which is, however, only of the order (M — m7r(TAf))~^/^, so that the total number 
of photons emitted is finite. In the case of a second order phase transition, the 
divergence disappears. However, since the temperature stays near Tm for finite 
time, the critical scattering peak remains visible. The threshold singularity of the 
cross section, which is responsible for the critical scattering effect, may be weakened 
due to higher order corrections in l/Nc- We nevertheless suspect that the cross 
section stays large at threshold. 

In the present work we have limited ourselves to the investigation of the photon 
pair production in the quark phase and did not take into account the contribution 
of the later expansion stage. This calculation is considered as a first step towards 
a more complete analysis of in medium effects on photon production in heavy ion 
collisions within chiral quark models. The hadronic tail of the evolution of a hot 
quark plasma phase also gives important contributions. A study of the processes in a 
pion gas, as 7r° — > 77 and tt+tt^ 77, is in progress p^. The applicability of quark 



models without confinement in the hadronic phase is limited to those kinematical 
situations, where unphysical quark production thresholds do not play a role. A first 
step towards the development of chiral confining quark models at finite temperature 
has been undertaken in Ref . . Also in these more realistic models the Mott effect 
occurs and we expect its manifestation in the two photon spectra. 
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APPENDIX A: MESON TO PHOTON DECAY VERTICES 

1. o" ^ 77 

The generic Feynman graph for the vertex ex — >■ 77 is given by Fig. ^ with F = 1. 
The contribution due to this diagram is 



- iKt^^ = -^T.J ^Tr [^Sf{q){-^Yqfe)^Sf{q - h^Sf^q + 

(Al) 

where we have introduced the symbolic four vector notation 

q = (iUn, q) h = (iu^, ki) /ca = h) . (A2) 

In Eq. ( |A1D , Un is a fermionic Matsubara frequency, whereas z/^ and ai are bosonic 
Matsubara frequencies belonging to the outgoing photons. It is understood that ium 
and iai have to be analytically continued to real values at the end of the calculation. 
The symbol Tr denotes the traces over spin, color and flavor degrees of freedom. 
Taking the trace over color and flavor yields a factor 



so that we have 



"-^^ 9 (3 (27r)3 [q2-nil][{q~h)^-ml][{q + k2y-ml] ' 

(A4) 
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where now the symbol tr denotes the spinor trace only. After taking the trace and 
noting that terms proportional to and fcg do not contribute to the final result due 
to the transversality of the photon, one has 



20 iV.mge^^ f d^q {k^k^ - gP^kik2) - gP\q^ - ml) + AqPq^ 



9 P 



(27r)3 [q2 - niMq - k.f - m2][(g + k,y - m^] 



(A5) 



Decomposing this integral and performing the analytical continuation leads to 

2ml\ 



{klyi-gf^k^k^) [l- 



Co{ki, -k2) 



97r '■^' "^ " V qiq2 ) 

where = 1/137 is the electromagnetic fine structure constant and 



(A6) 



Co{ki, -k2) 



167r2 



E 



d?q 



1 



(27r)3 [g2 _ rnf^ [{q - kiY - m^] [(g + ^2)^ - m? 



(A7) 



In the kinematical situation required throughout this work, ki 
\ki\ = \k2\ = this integral can be evaluated to be |19| 



-k2, k\ 



"-2 



r.^f^ , ptanh(f3E/2) , 
Co{ki, -k2) = -2V dp ^^^^ log 



E2(s -4^2) 



E-p 



E + p 



(A8) 



TT 



+ i-tanh(/3v/s/4) log 
s 



y/s - Js - Am? 



^+ Js- Ami 



and is thus a function of s only. Beside ([A6|) , there is a second contribution to the 
vertex from the crossed diagram of Fig. |^, which arises by exchanging ki ^ k2 and 
p ^ X. The contribution of this diagram is found to be equal to that of the first 
diagram. The function a{s) in Eq. ( p.l6|) may be thus written as 



ia{s) 



Qtt I 



-^]Co{k,,-k2 



(A9) 



This result is in agreement with the result of Ref. |28 . 



2. vr" 



77 



The vertex 7r° 77 can again be computed from Fi^ 
take r = i'y^T^. This leads to 



but now one has to 
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"^^^ P^J (27r)3 



(AlO) 



X 



Tr [iSf{q){-iYqfe)iSf{q - ki){i'^^T^)iSf{q + k2){-i'y^qft 



in the same notation as in Eq. (|AT|). The trace over color and flavor degrees of 
freedom yields now 



(All) 



and one has 



E 



d^q tr ii + mg)Yii-h + ^gh5{i + h + T^gW 



(A12) 

The only terms of the spinor trace which give a contribution are those, which involve 
75 and four other gamma matrices. One obtains 



(A13) 



and the same contribution from the crossed graph. In total, one has for b{s) required 
in Eq. (|2T7| ) 



(A14) 



This agrees with the result of Ref. 24 



APPENDIX B: SQUARED TRANSITION AMPLITUDES 



To compute the integrated cross section, one need the square of the transition 
amplitude, averaged over incoming and summed over outgoing states. 



1 



4Ar2 

c s,c 



J2 \Mt + Mu + Ms,a + Ms,n\ 



(Bl) 



We will quote our results for the single terms in Eq. (pl| ), beginning with the t- and 
M-channel diagrams. For the squares of these diagrams, one has 



Am 



Am 



s,c 



2 _32n^ 



2^4^,2 
em 



m{t - m2)2 
Nc{u — m^Y 



(m^ -u){s - 2m^) - {u + w?Y 



(B2) 
(B3) 
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whereas the interference term yields 
1 



Performing the angular integration, one obtains 



Cmin 



(B4) 



(B5) 



X 



(s + Am^s — 8m^) log 



s + x/s — Arrig 



y/s - Js - 4m2 



— (s + 4m — 4mg) 



and from this the result for the integrated cross section without resonances 
For the squares of the s-channel diagrams, one has 



4iV2 



a(s)P,(v^,0) 



(B6) 



and 



4iV2t^'- 4iV, 



6(s)P^(v^,0) 



(B7) 



The interference term of the scalar and pseudoscalar s-channels vanishes. For the 
interference terms of the s-channel amplitudes with the direct amplitudes, one ob- 
tains 



1 „ 27r(7facm^o 



-ia{s)V,{^s, 0) 



4Ar2 
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(B8) 



(B9) 



for the sigma exchange and 
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for the pion exchange. Integrating over t leads to the results 

1 9 ^' 



^min 



and 



+ s 



6(s)I)^(v^,0) 



dt—- {Ms,a + Ms,n) {Mt + Mu) 

^min TiiV^ 



(B12) 



(B13) 



2TTmqq'jacmS . 



log 



s — Am?q 



y/s - Js - 4m2 



From this, one immediately obtains the integrated cross section via Eq. (|2.18| ). 
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FIGURES 

FIG. 1. Mass spectrum of the SUf{2) NJL model. Solid line: pion mass, dashed line: 
double quark mass, dotted line: sigma mass. 

FIG. 2. Feynman diagrams for qq 77. Solid lines denote quarks, wavy lines photons 
and double lines mesons. Graph (a) shows the direct diagram, graph (b) the s-channel 
diagrams. Two more diagrams, which arise by the substitution k\ k2-, p \ are not 
shown. 

FIG. 3. Feynman diagram for the meson to photon decay vertex. Here one has to set 
F = 1 for cr ^ 77 and F = i^^T^ for tt^ — > 77. A second diagram arises by the substitution 
k\ <-> k2 and p <-> A. 

FIG. 4. Temperature evolution for the parameter set Tq = 250 MeV and tq = 1 fni/c. 
The solid line shows the first order scenario, the dashed and dot-dashed curves the second 
order scenario for AT = 10 MeV and AT = 1 MeV, respectively. The dotted line indicates 
the Mott temperature. 

FIG. 5. Comparison of photon spectra containing resonance contributions (solid line) 
and photon spectra without resonance contributions (dashed line) in the case of a first 
order transition at central rapidity as a function of the invariant pair mass M. 

FIG. 6. Comparison of photon spectra containing resonance contributions for a first 
order phase transition (solid line) and for a second order phase transition with various 
values of AT at central rapidity as a function of the invariant pair mass M. 
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Figure 2 



Figure 4 
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